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ABSTRACT 

Using stellar luminosity functions derived from ROSAT data, the contributions of Galactic 
stars to the diffuse X-ray background are calculated for ROSAT PSPC energy bands. The model 
follows that of Guillout et al. (1996), but uses i?0&4T' rather than Einstein data to determine the 
intrinsic luminosity distributions. The model adequately predicts the numbers of stellar sources 
observed in deep ROSAT surveys. The contribution of unresolved stellar sources to the ROSAT 
All-Sky Survey at the Galactic poles is 6.85, 4.76, and 4.91 x 10~ 6 counts s _1 arcmin -2 in bands 
R12 (j keV), R45(| keV), and R67(1.5 keV), respectively, which is equivalent to 4.66, 31.3 and 
26.9 xlO -14 ergs cm~ 2 s" 1 deg~ 2 . 

Subject headings: X-rays-diffuse background — X-rays:stars — X-rays:ISM 



1. Introduction 



1.1. Motivation 



Recent advances in the study of the soft X- 
ray background have attempted to place limits 
on the amount of emission from diffuse extra- 
galactic sources. Cen & Ostriker (1999) suggested 
that such a background might exist in the \ keV 
and | keV range as a result of gas accretion on 
sub-cluster size objects forming filamentary struc- 
tures on cosmological scales, the exact form of 
these structures being sensitive to the cosmological 
model. At higher energies (above ~ 1 keV) truly 
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diffuse cxtragalactic X-ray emission is ruled out as 
such gas would produce unacceptably large per- 
turbations in the microwave background (Wright 
et al. 1994). Such considerations do not hold at 
2 keV, but the multiplicity of backgrounds (at 
least two diffuse Galactic backgrounds, as well as 
unresolved Galactic point sources, and unresolved 
extragalactic sources) make the detection of a dif- 
fuse extragalactic background difficult. ROSAT 
observations have been used to attack the prob- 
lem on several different fronts. 

The extragalactic background due to unre- 
solved point sources can be characterized as a 
power-law at E>1 keV, and a variety of results, in- 
cluding recent Chandra observations (Mushotzky 
et al. 2000) suggests that this power-law extends, 
without a break, to lower energies. Snowden et al. 
(2000) and Kuntz & Snowden (2000), following 
the work of Snowden et al. (1998), have stud- 
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ied the emission that is not due to unresolved 
extragalactic sources. Snowden et al. (1998) sepa- 
rated the \ keV emission originating from beyond 
the Galactic absorption from that due to the Lo- 
cal Hot Bubble, and this work was extended by 
Kuntz & Snowden (2000) to the remainder of the 
ROSATband. Although Kuntz & Snowden (2000) 
characterized the emission originating beyond the 
Galactic absorption in the f keV band, they were 
not able to effect a separation of the emission 
from beyond the Galactic absorption into its var- 
ious components: the Galactic halo, the diffuse 
extragalactic emission, and unresolved Galactic 
point sources. Thus, although there is almost no 
doubt that there is a Galactic halo, the upper 
limit to the diffuse extragalactic emission in the 
| keV band remains the total emission originating 
beyond the Galactic absorption. Included in the 
measure of the emission originating beyond the 
Galactic absorption was the bulk of the contri- 
bution due to unresolved Galactic point sources. 
It is thus necessary to calculate that contribu- 
tion to further tighten the limits on the diffuse 
extragalactic background. 

From the model of Cen & Ostriker (1999), 
the diffuse extragalactic background will not be 
smooth, and the dumpiness is a function of the 
cosmological model. Kuntz (2000) sought to place 
limits on the spatial fluctuations that might be 
due to the diffuse extragalactic background, thus 
it was necessary to determine the degree of fluctu- 
ation due to unresolved Galactic point sources. 

It is thus particularly interesting for the study 
of the diffuse extragalactic background to study 
the luminosity function of Galactic point sources, 
which, at high Galactic latitudes, are stars. 

1.2. Previous Work 

A stellar log N-log S relation in the form of a 
cumulative luminosity function (CLF, the number 
of sources, N, with fluxes greater than S) has been 
calculated by Guillout et al. (1996) from a mixture 
of Einstein IPC and ROSATPSPC data, using the 
Bienayme et al. (1987) model of the spatial distri- 
bution of stars. The Guillout et al. (1996) study 
concentrated on comparing the model distribution 
of stellar sources with observed count rates as a 
function of position, the complement of the infor- 
mation required for studies of the diffuse X-ray 
background. 



There are a number of difficulties in applying 
the published Guillout et al. (1996) CLF to the 
problem of unresolved stellar sources. First, the 
Guillout et al. (1996) study published the CLF 
for the entire ROSATPSPC band-pass. Since the 
angular structure of the diffuse X-ray background 
is strongly energy dependent, studies of the dif- 
fuse X-ray background are generally confined to 
only some portion of that band-pass. Given the 
strong differences in the absorption by the ISM 
from one end of the band-pass to the other, no 
simple scaling can be applied to the published re- 
lations to obtain the relation for any sub-band. 
Second, as can be seen in Guillout et al. (1996), 
and again in our Figure 1, although older stars 
have lower X-ray luminosities, they form such a 
large fraction of the total stellar population that 
they dominate the luminosity function, especially 
at lower flux levels. However, this is the popu- 
lation that is the least well characterized in the 
Guillout et al. (1996) study; they were forced to 
take data from Einstein IPC studies of field stars. 
Far superior ROSAT studies have since become 
available (Schmitt et al. 1995; Schmitt 1997) 

The following model calculation uses the best 
current data to derive the luminosity function at 
high Galactic latitude, which will be of greatest 
interest to those studying the diffuse extragalactic 
background. Rather than the CLF, we calculate 
the more useful form, 

N(S)d\ogS, (1) 

which allows the calculation of 

f S SN(S)d\ogS, (2) 
Jo 

the flux in unresolved sources. The model is, in 
principle, capable of calculating the stellar contri- 
bution at low Galactic latitudes, but some poorly 
determined input parameters, the intrinsic lumi- 
nosity probability distribution for young stars, 
greatly increase the uncertainties near the Galac- 
tic plane. The irregularity of the spatial distribu- 
tion of the ISM also introduces some uncertainty, 
more in softer bands than in harder, and more at 
lower Galactic latitudes. 

2. Data 

In order to calculate the contribution of stel- 
lar sources, the stellar populations were binned 
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Fig. 1. — The luminosity function, N(S)d\ogS, for band R17 shown with the fractional contributions of the 
different age classes (left) and spectral types (right). 



by spectral type St and age A. We used the 
spectral-type divisions of Schmitt et al. (1995) and 
Schmitt (1997) and the age divisions of Guillout 
ct al. (1996), all of which arc listed in Table 1. The 
division by age is as follows: the 0-0.15 Gyr class is 
represented by the Pleiades (171 ROSAT detected 
members, ~ 0.1 Gyr, 127 pc), the 0.15-1.0 Gyr 
class is represented by the Hyades (191 detected 
members, <~ 0.6 Gyr, 45 pc), and the 1.0-10.0 Gyr 
class is represented by the field population. 

The data arc insufficient to determine the X- 
ray luminosity functions for subsamples of the 1-10 
Gyr field population, which, especially for the later 
type stars, will be nearly uniformly distributed 
through the age range. However, stars in the 1-2 
Gyr age range will have a very different z-height 
distribution than stars in the 5-10 Gyr age range. 
Thus, for calculation, the stars in the 1-10 Gyr age 
class are subdivided into the age ranges used by 
Bienayme ct al. (1987) (see Table 2) although the 
same X-ray luminosity function is used for all the 
age sub-divisions in the 1-10 Gyr age class. 

A number of other well-studied clusters could 
be included with the Hyades for a sample of stars 
in the 0.15-1.0 Gyr group; Coma Berenices (42 de- 



tected members, ~ 0.5 Gyr, 80 pc) and Praesepe 
(56 detected members, ~ 0.6 Gyr, 160 pc). X-ray 
emission from Coma Berenices members is simi- 
lar to that from Hyades members (Randich et al. 
1996), while that from Praesepe members is lower 
than that from (coeval) Hyades members (Randich 
& Schmitt 1995). Given that the Hyades would 
seem to be representative of its age class, and that 
there is a relatively larger sample available from 
the Hyades cluster than from other clusters, the 
other clusters have not been included in this age 
class. However, it should be borne in mind that 
the Hyades luminosities may need to be adjusted 
downwards to truly represent its age class. 

Although X-ray emission has been detected 
from nearly every spectral type, the dominant con- 
tribution by stellar sources is from low-mass main 
sequence dwarfs, due both to intrinsic luminosity 
and the number of sources. 

Stars of spectral type A have been excluded 
from the model. The intrinsic X-ray luminosity 
probability distribution is poorly determined for 
the younger age classes, and is currently unde- 
termined for the 1-10 Gyr age class, which would 
dominate the emission from the A spectral type. 
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The intrinsic X-ray luminosities, where measur- 
able, are about an order of magnitude below those 
of F stars, and the mid-plane number density is 
lower than F stars, by a factor of 0.005 for the 
1-10 Gyr age class. Since the younger age classes 
will have low z-heights, the error introduced will 
be minimal for high Galactic latitudes. Guillout 
ct al. (1996) estimate the maximum contribution 
from A stars to be ~ 2%. 

White dwarfs have also been excluded from the 
calculation as their contributions are typically two 
orders of magnitude below that of main-sequence 
stars in the \ keV band, and significantly lower in 
the | keV and 1.5 kcV bands. 

RS CVn systems have also been excluded from 
the calculation. Although they have high lumi- 
nosities, the median intrinsic luminosity, log L x ~ 
30.3 (here, as in the remainder of this paper, L x 
is in ergs s^ 1 in the 0.1-2.4 keV band) (Dempsey 
et al. 1993), their number density is quite low, 
probably less than 1 x 10~ 4 pc~ 3 (Ottmann & 
Schmitt 1992) as is their scale height (~ 325 pc). 
At high Galactic latitudes, any RS CVn system 
is likely to be close enough to be detected and 
not contribute to the "diffuse" background. Gi- 
ants and super giants have been excluded from 
the calculation as their proper calculation would 
require extensive analysis beyond the scope of this 
paper. Nearly all giants earlier than K3 are X-ray 
sources with log Lx ~ 27.3, similar to solar type 
stars (Schroder et al. 1998; Hiinsch et al. 1998) but 
the source population is much smaller, and more 
closely restricted to the plane. 

We have also excluded Population II stars from 
the calculation. If a comparison between disk RS 
CVn systems and Population II binaries is any in- 
dication of the relative emission from single Pop- 
ulation I and Population II stars, then the Popu- 
lation II stars should have emission that is more 
than an order of magnitude below that of Popula- 
tion I stars (Ottmann et al. 1997). 

2.1. Model Inputs 

For each spectral- type/age-class bin, the contri- 
bution to the luminosity function is determined by 
four parameters: fj, and ct, the mean and dispersion 
of the intrinsic luminosity probability function, c, 
which describes the z distribution of the stars, and 
no, the mid-plane (z = 0) density of stars. 



Following Schmitt et al. (1995) and Schmitt 
(1997), we have assumed that the intrinsic lumi- 
nosity function for a spectral- type/age-class bin is 
given by the log-normal distribution, 



P(logL x ) = 



— == exp — - 
ctv/21F 2 



\ogL x - n 



(3) 



The parameters [i(St,A) and (t(St,A) were de- 
rived from the i?OS^4T data sets listed in Table 1. 

Since these data sets contain points with upper 
limits, the parameters were determined by fitting 
the Kaplan-Meier estimator (Schmitt 1985) of the 
cumulative distribution function, where the uncer- 
tainty in the Kaplan-Meier estimator was taken 
from Fiegelson & Nelson (1985). The uncertain- 
ties were derived using a bootstrap method (Babu 
& Feigelson 1996). We corrected for binaries in the 
same manner as Schmitt et al. (1995): we have di- 
vided the flux equally between members of binary 
systems. Guillout et al. (1996) showed that cor- 
rection for binaries made little difference to the 
luminosity function. 

Since there were no late M stars detected in 
the surveys of the Hyades or Pleiades, we have 
extrapolated the u for the the 0-0.15 Gyr and 0.15- 
1.0 Gyr age classes from the fi for the 1-10 Gyr 
age class using the general trend observed in the 
other spectral classes; a Alogix = 0.92 decrease 
between the 0-0.15 Gyr and 0.15-1.0 Gyr classes, 
and a AlogLx = 0.98 decrease between the 0.15- 
1.0 Gyr and 1.0-10.0 Gyr classes. For the late M 
stars the same value of ct was used for each age 
class, that of the oldest age class; this tends to 
increase the number of sources with high fluxes, 
but not significantly. 

To convert from the listed X-ray luminosities 
to count-rates, we used the inverse of the con- 
version used by the original authors. Stern et al. 
(1995) used the conversion 6.0 x 10~ 12 ergs cm~ 2 



Stauffer et al. (1994) used 

2 „-l 



8.4 x 10 12 ergs cm * s 



s _1 [counts s _1 
a conversion of 
[counts s -1 ] -1 (for Nh = 0), characteristic of 
logT = 7.06. Due to the shape of the conver- 
sion curve, this is also the conversion factor for 
logT = 6.88. Schmitt et al. (1995) used the con- 
version factor 4 (5.30.072 + 8.31) x 10" 12 ergs cm~ 2 



4 This relation was derived for a Raymond & Smith model 
two component corona fitted to a small number of high S/N 
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Table 1 
Stellar Data 



Type 


B-V 




No. a 


0-0.15 Gyr 
„b 


a b 


No. a 


0.15-1 Gyr 
M b 


<7 b 


No. a 


1-10 Gyr 




A d 


-0.01-0.29 




16/39 


28.51±0.26 


0.99±0.16 


10/35 


27.84±0.39 


0.93±0.29 


16/62 6 






F 


0.3-0.6 




30/33 


29.39±0.07 


0.32±0.06 


57/60 


28.93±0.05 


0.28±0.03 


28/29 


28.00±0.14 


0.96±0.14 


G 


0.61-0.8 




19/26 


29.39±0.20 


0.67±0.18 


27/33 


28.64±0.12 


0.39±0.05 


29/38 


27.36±0.17 


0.82±0.13 


K 


> 0.8 


< 8 


40/44 


29.46±0.07 


0.42±0.04 


51/86 


28.29±0.06 


0.45±0.08 


17/17 


27.62±0.17 


0.48±0.12 


M (early) 


> 0.8 


8-15 


44/55 


29.19±0.05 


0.32±0.04 


56/174 


27.89±0.14 


0.72±0.12 


66/78 


26.86±0.10 


0.77±0.10 


M(latc) 


> 0.8 


> 15 





29.12 f 


0.49 f 





28.20 f 


0.49 f 


16/19 


27.22±0.22 


0.49±0.24 



Spatial Distribution 


Bienaymc ct al. Eq. Al 


Bienaymc ct al. Eq. A2 


Bienaymc ct al. Eq. A2 


Stellar Sample 


Pleiades 


Hyades 


Field Stars 

d<13 pc for F & G stars 
d<7 pc for K & M stars 


Data 8 


Stauffer et al. (1994), PSPC 


Stern et al. (1995), RASS 
Pyc et al. (1994), PSPC 


Schmitt ct al. (1995), PSPC 
Schmitt (1997), PSPC 



a The number of detected stars in the sample divided by the total number of stars in the sample. 

b The quantities /i and a are in log 10 L x where L x is in units of ergs s _1 in the 0.1-2.4 keV band. The uncertainties were determined as 
described in the text. 

c Sincc Schmitt (1997) and Schmitt ct al. (1995) convert from PSPC counts to luminosity using (5.30HR + 8.31)xl0~ 12 erg cm~ 2 s" 1 [count 
s -1 ] -1 where HR is a hardness ratio, and since HR is luminosity dependent, the a derived from Schmitt (1997) and Schmitt et al. (1995) data 
are likely to be systematically larger than those derived from other data sets where the PSPC count to luminosity conversion was made with a 
standard 6 X 10~ 12 erg cm~ 2 s _1 [count s -1 ] -1 . This difference in conversion may also introduce a very small offset between these luminosities 
and those of the other age classes. Our fit values were similar to those of Schmitt ct al. (1995) and Schmitt (1997), and we have retained their 
values. 

d Spectral- type A is not used in the model calculations. Note the large uncertainties in the parameters of the intrinsic luminosity distribution. 
The parameters could not be calculated for the field sample as the detected fraction for a sample of main-sequence A stars within 30 pc is ~ 1/3, 
and upper limits for the non-detections have not been published. 

"Number in the RASS Point Source Catalogue (Hiinsch et al. 1999). In the r < 13 pc sample of Schmitt (1997) there are four A stars, none 
of which were detected. 

f Extrapolated from the 1.0-10.0 Gyr class due to lack of data. See text. 

g "PSPC" refers to data obtained from pointed PSPC observations, "RASS" refers to data extracted from the ROSAT All-Sky Survey. 



Table 2 
Stellar Data 



Type Age Group (Gyr) 





0-0.15 




0.15-1 




1-2 




2-3 




3-5 




5-10 




A 


2.16 x 10" 


-4 


3.93 x 10" 


-4 


1.01 x 10" 


-4 


2.55 x 10" 


-7 












F 


2.93 x 10 


-4 


6.14 x 10" 


-4 


4.40 x 10" 


-4 


2.59 x 10" 


-4 


2.29 x 10" 


-4 


6.99 x 10" 


-5 


G 


2.78 x 10" 


-4 


5.82 x 10" 


-4 


4.17 x 10" 


-4 


2.88 x 10" 


-4 


4.39 x 10" 


-4 


9.15 x 10" 


-4 


K 


3.31 x 10 


-4 


6.93 x 10" 


-4 


4.96 x 10" 


-4 


3.42 x 10" 


-4 


5.22 x 10" 


-4 


1.18 x 10" 


-3 


M (early) 


4.28 x 10" 


-3 


8.96 x 10" 


-3 


6.42 x 10" 


-3 


4.43 x 10" 


-3 


6.75 x 10" 


-3 


1.53 x 10" 


-2 


M(late) 


2.45 x 10 


-3 


5.12 x 10" 


-3 


3.67 x 10" 


-3 


2.53 x 10" 


-3 


3.86 x 10" 


-3 


8.73 x 10" 


-3 


c a 


0.0140 




0.0279 




0.0457 




0.0662 




0.0867 




0.0958 





a Paramcter describing the cllipticity of the stellar distribution. See equations Al and A2 of Bicnaymc 
ct al. (1987). 



s 1 [counts s : ] 1 , where HR is denned by (R4 7- 
R82)/(R47+R82) and RXX is the count rate in 
band XX and the bands are defined in Snowden 
et al. (1997) and again in Table 3. 

From the band R17 count rates it is neccessary 
to derive count rates in narrower bands for this 
study to be of use. To divide the band R17 into the 
R12 and R47 bands, one can use the log L X -HR 
developed by Schmitt (1998). Schmitt suggested 
that for each spectral type, HR is a linear function 
of the log of the luminosity, 

\ogL x =aHR + b; (4) 

the slopes of these relations seem to be the same 
for all spectral types, but the intersections are dif- 
ferent. The parameters of these relations were 
determined for the field stars from the data of 
Schmitt et al. (1995) and Schmitt (1997) by Kuntz 
(2000). 

The applicability of these relations to younger 
stars is not clear. Extracting the Hyades members 
identified by Stern et al. (1995) from the ROSAT 
Bright and Faint Source Catalogues (Voges et al. 
1999, 1996), we found that the field star relations 
fit the Hyades stars well for spectral types F and 
G. For the M stars, the slopes of the relations ap- 
peared consistent, but the field star offsets were 
too large (i.e., too hard) for the Hyades stars. 
Since there is a substantial dispersion around the 
relation for both the field stars and the Hyades, 

spectra. See Fleming et al. (1995) for a complete descrip- 
tion. 



and there may be significant biases in the selec- 
tion of sources for the Bright and Faint Source 
Catalogues, we have used the field star relations 
for the younger stars as well. There are not enough 
Pleiades stars in the Faint Source Catalogue to at- 
tempt a similar analysis. 

To further subdivide the flux, we investigated 
the relation between HR and the spectral shape. 
From the HEASARC archive of ROSAT PSPC 
exposures we extracted the spectra of all single 
field dwarfs of spectral classes A through M that 
had more than <~ 400 counts, about 16 stars. 
To these spectra we used XSPEC (Arnaud 1996) 
to fit two Raymond & Smith (1977) components 
with temperatures Ts and Th and normalizations 
N s and N H . We found that the temperature of 
the softer component, Ts, did not vary with HR 
(log Ts ~ 6.23), but that for stars in an individual 
spectral type, the temperature of the harder com- 
ponent, Th, varied linearly with HR 5 . Given the 
Th for a given HR, we then solved for the Ns/Nh 
ratio that would produce the correct value of HR. 
This relation between Th and HR restricts the 
range of allowable HR to > —0.86 and < 0.37 (for 
M stars) and < 0.22 (for F and G stars). 

We checked this T H , N S /N H , HR relation 

5 We found that the F and G stars followed one T H -HR 
relation while the M stars followed another. There were 
only two K stars in the sample, one of which fell on the 
relation for F and G stars, the other of which fell between 
the FG and M relations. Of the M stars, half were dM and 
half were dMe, so the difference between the FG and M 
relations is not due to selecting the brighter M stars. 
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Table 3 
Broad Band Definitions 



Band 


PI Channel 


Encrgy a 
(keV) 






Rl 


8-19 


0.11-0.284 


1 R12 




R1L = R8 


11 - 19 


0.11 - 0.284 


> or 


= \ keV 


R2 


20 - 41 


0.14-0.284 


1 R1L2 




R3 


42 - 51 


0.20-0.83 






R4 


52 - 69 


0.44-1.01 ] 


R45 


= ! keV 


R5 


70 - 90 


0.56-1.21 1 


R6 
R7 


91 - 131 
132 - 201 


0.73- 1.56 1 
1.05-2.04 J 


> R67 


= 1.5 keV 



Note. — The notation band RXY denotes the sum of all of the bands 
between (and including) RX and RY. In the ROSAT source catalogues, 
HR1=(R47-R12)/(R47+R12) and HR2=(R67-R45)/(R67+R45). 

a 10% of peak response. 



against the Hyades stars extracted from the 
Bright and Faint Source Catalogues using (R47- 
R82)/(R47+R82) vs. (R67-R45)/(R67+R45) 
graphs, and found that within the large uncer- 
tainties of the data, the relations derived from 
field stars also described the Hyades stars. 

Thus, for a given X-ray luminosity, we used the 
Schmitt relation between the X-ray luminosity and 
the hardness ratio to determine HR, and then use 
these relations between Th, Ns/Nh, and HR to 
derive a spectral energy distribution that can then 
be used to divide the total number of counts into 
the counts for each of the ROSAT bands. Since 
the log L X -HR relation can produce values of HR 
ouside the range allowed by the T H , N S /N H , HR 
relation, values of HR outside this range were re- 
placed by the closest extreme. 

The spatial distribution of each spectral- 
type/age-class was assumed to be given by equa- 
tions Al and A2 of Bienayme et al. (1987), where 
c(A), the axial ratio for ellipsoids describing the 
disk populations, were taken from the "best fit" 
column of their Table 1 and are listed in our Ta- 
ble 2. 

To maintain some consistency between the ini- 
tial mass function (IMF), the stellar formation 
rate (SFR), and the cr z -Age relation, the mid- 
plane density of each spectral- type/age-class bin, 
no (St, -A), was derived in a manner similar to that 



used by Guillout ct al. (1996). The disk was as- 
sumed to have an age of 10 Gyr, a constant SFR, 
and an invariant IMF. The IMF was taken from 
Haywood et al. (1997): 



where 



dN = m-^+^dm (5) 



x = 0.7 for m < m & 

x = 1.5 for Itoq < m < 3m (6) 

x = 2.0 for m > 3m . 

The total number of stars was normalized by com- 
paring the number produced by the above formula 
with r < 7 pc in the 8 < My < 15 range (i.e., 
the K and early M stars) with the number actu- 
ally observed (79). The fraction of the number 
of stars in each age-class was calculated using the 
main-sequence lifetimes listed in Scalo (1986) and 
the mass-absolute magnitude-spectral type rela- 
tions of Kroupa et al. (1993). 

2.2. Method 

For each spectral- type/age-class bin, the in- 
trinisic luminosity distribution was evaluated for 
AlogLx = 0.1 intervals. For each distance in- 
terval, Ar = 10 pc, from to 3000 pc, the 
value of log Lx was adjusted for r~ 2 diminution, 
absorption by the ISM, and converted to count 
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Fig. 2. — The relative number densities for each 
age group as a function of z- height, compared to 
the absorbing column through which they are ob- 
served. Top: The total absorbing column through 
which and object of z-height R is observed. The 
t = 1 point for band R12 is marked; to achieve the 
same t for band R47 requires a much greater col- 
umn density. Bottom: The relative number den- 
sities for stars of a given age group at a z-height 
R. Each curve has a mid-plane number density of 
unity. 

rates in a given band. The probability of finding 
a star with flux S in a given band, as a func- 
tion of S, was then interpolated to a standard 
A log (S'/counts s' 1 ) = 0.1 spacing, and scaled by 
the size of the volume clement and the relative 



stellar density for that z-height. The luminos- 
ity function for each spectral- type/age-class was 
then summed over z-height, and scaled by the to- 
tal number of stars 

The vertical distribution of the ISM was taken 
from Lockman (1984), and was assumed to be 
plane-parallel. The absorption cross-sections were 
taken from Morrison & McCammon (1983). The 
absorption was calculated using the spectral en- 
ergy distribution described above. Although the 
bulk of the stars lie beyond the bulk of the absorp- 
tion (see Figure 2), only the R12 band reaches 
column densities for which r = 1. (The next 
most sensitive band is R17, for which r = 1 at 
~ 5.5 x 10 20 cm -2 .) Thus, the luminosity func- 
tion is not too sensitive to absorption by the ISM 
at high Galactic latitudes. 

The aggregate luminosity function for band R17 
is shown in Figure 1 for the north Galactic pole. 
As might be expected, the luminosity function is 
dominated by the older age classes because, al- 
though their intrinsic luminosities may be an or- 
der of magnitude smaller than the younger age 
classes, they account for the bulk of the number 
of sources. For similar reasons, the luminosity 
function is dominated by the longer-lived spectral 
classes. The high-count-rate end of the luminos- 
ity function is strongly influenced by the younger 
age-classes as they have larger intrinsic luminosi- 
ties and smaller scale-heights. 

Figure 3 shows the effect of the uncertainties 
of the parameters on the luminosity function (top 
third), unresolved flux (middle third), and cumu- 
lative luminosity function (bottom third) for a sin- 
gle spectral- type/age-class (early M, 1-2 Gyr) in 
band R47 (the most interesting band for cosmolog- 
ical purposes). In this band, the effects of Galac- 
tic absorption are minor. The parameter uncer- 
tainties are their typical formal uncertainties. For 
typical survey limits, -4 < \ogS(R47) < -2, the 
predicted unresolved flux is uncertain to a factor 
of two due to the formal uncertainties. The largest 
sources of uncertainty are the \x and a of the in- 
trinsic luminosity distribution function. Although 
the b paramter of the log L X -HR relation also ap- 
pears to contribute a significant uncertainty, the 
true uncertainty is likely smaller than the value 
shown here, which includes what may be a signif- 
icant intrinsic scatter around the relation. Uncer- 
tainties in c cause little change. 
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Fig. 3. — The effect of the formal uncertainties of 
the input parameters on the luminosity function 
(top), the unresolved flux (middle), and the cu- 
mulative luminosity function (bottom), shown as 
the ratio between the luminosity functions. For 
the M(early), 1-2 Gyr age class we show the R47 
band luminosity function calculated with a sin- 
gle parameter increased by A, where the A is the 
typical formal uncertainty, divided by the fiducial 
luminosity function. 



The true uncertainty in the amount of unre- 
solved flux is much less than that shown in Fig- 
ure 3; as will be seen, the model can be checked 
against the observed N(S) over several orders of 
magnitude in source flux, and good agreement be- 
tween the model and the observed N(S) implies a 
smaller error in the unresolved flux. 

3. Comparisons 

The current model agrees fairly well with that 
of Guillout et al. (1996) except for the "late M" 
class. For this class we find a factor of <~ 10 more 
sources. Only a small amount of this difference can 
be attributed to the difference in luminosity dis- 
tributions; the luminosity distribution used here 
(/i = 27.22,cr = 0.49) is surprisingly similar to that 
derived from Einstein data (fi = 27.7, o ~ 0.42) by 
Barbcra ct al. (1993). It is not clear whether the 
difference in number counts for the late M class is 
due to a different IMF, a different low-mass cut- 
off for X-ray emission (0.08 M was used here), 
or merely a different mass-to-spectral-type conver- 
sion. All of these quantities are, of course, quite 
uncertain for the low mass end of the main se- 
quence. 

We have checked our model against the num- 
ber of stellar point sources found in surveys and 
against the stellar log N-log S relation of Krautter 
et al. (1999). These comparisons are listed in Ta- 
ble 4 and shown in Figure 4. The predicted num- 
ber listed in the Table is the number of sources 
with flux above the survey limit at the (I, b) posi- 
tion of the survey, in degrees -2 , multiplied by the 
survey area and the fraction of sources for which 
identifications have been made (the "ID complete- 
ness" ) . 

As can be seen, the model predicts well the 
band R47 source counts for fluxes below log 10 S ~ 
—2.5. At higher fluxes we can compare our model 
to the wide-field (each field has ~ 144 square de- 
grees) data of Krautter et al. (1999) which is com- 
plete in band R17 to 0.03 counts s _1 , (and thus 
must be complete to at least that level in the sub- 
band R47). We find that in band R47, that for 
their field IVac (b = 79°), the model is consis- 
tent with observations for fluxes below 0.1 counts 
s _1 . However, below that threshold, the model 
over-predicts at low latitudes (30° < \b\ < 40°, 
their fields I, II, III, and VI). Above this thresh- 



Table 4 
Predictions 



Survey 


(t,b) 


Band 


Area 


ID 


Source 


Number 


Number 




2000 




dcg 2 


Completeness 


Limit 


Predicted 


Found 












10~ 3 c/s 






Bower et al. (1996) 


( 96.4, 29.8) 


47 


0.21 


0.90 


.763 


3.07 


3 


Schmidt ct al. (1998) 


(149.5, 53.5) 


47 


0.28 


1.00 


.459 


3.71 


3 


Zamorani ct al. (1999) 


(270.2,-51.8) 


47 


0.20 


0.84 


.309 


3.73 


4 


McHardy et al. (1998) 


( 85.0, 75.9) 


47 


0.17 


1.00 


.175 


3.84 


5 


Carballo et al. (1995) 


(162.6, 51.1) 


47 


2.20 


0.57 


1.42 


11.11 


9 


Boyle et al. (1997) 


various 


47 


3.92 


0.89 


1.7 


18.93 


20 


Comparison to measured log N-log S function 


Krautter et al. (1999) 


various 


17 
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0.6 a 


0.64 
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170 


0.058 


0.057 
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1700 


0.0032 


0.0030 



a Value from the log N-log S function graphed in their paper. We have assumed that this function is valid for 
b = 90°. 
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Fig. 4. — The predicted cumulative luminosity 
function with measured values, x: values in band 
R47 from the deep surveys listed in Table 4, +: 
(shown with error bars) values in band R47 from 
the field IVac (b = 79°) source list of Appcnzcllcr 
ct al. (1998). □: (not shown with error bars) 
values in band R47 from fields I, II, III, and VI 
of Appenzeller et al. (1998). These fields have 
30° < |fe| < 40°. See text for further explanation. 



old, the model under-predicts the source counts at 
all latitudes. A similar trend is seen in the band 
R17 source counts; the model somewhat under- 
predicts the counts at high Galactic latitudes and 
over-predicts at low Galactic latitudes for fluxes 
greater than 0.03 counts s _1 . Given the shape of 
the luminosity functions (Figure 1), this error is 
likely due to the poorly determined \i and a of the 
youngest age group. (Nota Bene: the caveat given 
in §2 concerning the Hyades luminosity function, 
as well as our use of a log L X -HR relation for the M 
stars that may be too hard, decreasing the effect 
of absorption, §2.1.) 

The problems with the model at bright fluxes 
do not cast serious doubt upon our prediction at 
lower fluxes; as can be seen from Figure 3, the 
brighter fluxes are much more sensitive to uncer- 
tainties in fi and a, and the brighter fluxes are 
much more sensitive to the younger age groups. It 
should also be noted that the uncertainties in the 
unresolved flux at the brighter fluxes will be sig- 
nificantly smaller than those of the source counts. 

4. Model Description 

It is instructive to see how the luminosity func- 
tion changes with detection band, Galactic lati- 
tude, and Galactic longitude. 

Figure 5 shows the stellar luminosity function 
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Fig. 5.— The luminosity function for different 
bands in the direction of the north Galactic pole. 
Solid line: band R17, +: band R12, □: band R45, 
O: band R67, Dotted line: band R47 (almost in- 
distinguishable from band R45). 

at the north Galactic pole for several ROSAT 
bands. The R12 band is substantially different 
from the R45 band due both to the different in- 
strumental response (which produces a horizontal 
shift in the curves) and the differing effects of ab- 
sorption (which increases the number counts be- 
low the peak and decreases them above the peak, 
though not at the highest fluxes). The anomolous 
behavior of band R67 is due to the discontinuity of 
the log L X -HR relation at HR = —0.86; the rate of 
decrease of the band R67 flux with decreasing HR 
is greater than for the other bands, thus the dis- 
continuity is more pronounced for band R67 than 
for the other bands. As the anomolous behavior 
occurs two orders of magnitude below the peak of 
the luminosity function, it does not cause signifi- 
cant errors. 

Figure 6 shows the ratio of the stellar luminos- 
ity function at b = 30° to that at the Galactic pole 
for the R12 and R47 bands. We have assumed a 
plane-parallel absorption model. For typical deep 
ROSAT surveys, there is a factor of ~ 2.5 greater 
number of sources at b = 30° than at b = 90° 
in band R47. The increase in number counts for 
band R12 is less severe due to absorption. As can 



-6 -4 -2 

Log, (S) 

Fig. 7. — The luminosity function for main- 
sequence stars, white dwarfs, and extragalactic ob- 
jects, for band R12 (dotted) and band R47 (solid). 
The relations are shown for a total Nh — 5 x 10 19 
cm~ 2 (further to the right) and for Njj = 2.1 x 
10 20 (further to the left). For band R47, the re- 
lations for the different values of Nh are indistin- 
guishable. The band R47 relation for white dwarfs 
is not shown, but would be significantly below the 
R12 relation. The point source detection limit for 
the R OS A T All-Sky Survey is also shown. 

be seen from the Figure, at b = 30° the effect of 
the Galactic radial density gradient is negligible. 

5. Results 

Figure 7 shows the luminosity functions in 
PSPC bands R12 (dotted) and R47 (solid) for 
main-sequence stars, white dwarfs 6 , and the extra- 
galactic luminosity function derived by Hasingcr 
ct al. (1998). For band R47, above the RASS 
detection limits, stellar sources account for one 
third to one half of point sources, in agreement 
with studies such as those listed in Table 4, while 
for band R12, the fraction depends upon the total 



The white dwarf luminosity function for single white dwarfs 
is shown to demonstrate its inconsequcntiality. It was cal- 
culated using the data of Fleming et al. (1996) with the 
assumption that the intrinsic X-ray luminosity of white 
dwarfs extends down to at least logLx = 27.5 
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Fig. 6. — The ratio of the luminosity function at b = 30° to the luminosity function at b = 90°. Solid: The 
ratio for £ = 90°, Dotted: the ratio for I = 0°, Dashed: the ratio for £ = 180°. 



absorbing column. 

Figure 8 shows the contribution to the RASS of 
unresolved stellar sources as a function of Galac- 
tic latitude. Table 5 shows the model prediction 
for the contribution of the unresolved Galactic 
stars to the diffuse X-ray background at the north 
Galactic pole for N H = 2.1 x 10 20 cm" 2 . The con- 
version of count rates to energy use the spectral 
shape of the unresolved flux and the energy ranges 
listed in the table. At the RASS point source 
detection limit in band R45 the unresolved flux 
has R12:R45:R67 =1.51:1:1.09, and those ratios 
are similar for the unresolved flux at the RASS 
point source detection limits in the other bands. 
The unresolved flux in the R47 band (0.44-2.04 
kcV) is 3.64(3.51) x 10~ 13 erg cnT 2 s" 1 dcg~ 2 at 
the RASS band R45(R67) point source detection 
limit. 

In bands R12, R45, and R67, the unresolved 
stellar sources can account for, at most, 2%, 10% 
and 51% of the flux that is thought not to be due 
to the Local Hot Bubble or the unresolved extra- 
galactic point sources (essentially the component 
that Kuntz & Snowdcn (2000) attribute to the 
Galactic halo, but which can be distributed be- 



tween the Galactic halo and diffuse extragalactic 
emission). The variation of the contribution from 
unresolved stellar sources with Galactic latitude 
would be unnoticeable in the RASS given the un- 
certainty in the RASS and, more importantly, the 
variable effects of Galactic absorption. 

The deep band R47 surveys discussed in §3.2 
suggest that the predicted unresolved flux for 
surveys at high Galactic latitudes with limits 
logS*(i?47) < —2.5 is correct. At higher fluxes 
the comparison of the model with observations 
is more equivocal, but the high Galactic latitude 
source counts in band R47 do seem to be correct, 
(at least for S < 0.3 counts s _1 ) and thus the 
unresolved flux should be correct. Due to the un- 
certainties introduced by the younger stars, the 
model is not adequate at lower Galactic latitudes 
(|6| < 30°), particularly for the high flux source 
counts as these are dominated by younger stars. It 
should also be noted that both the younger stars 
and the Galactic absorption will be distributed in 
a much more clumpy manner than used by this 
model, so a great deal of field to field variation is 
to be expected. 
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Table 5 
Unresolved Stellar Flux 



ROSAT Energy RASS Unresolved Flux 
Band Rangc a PSDL b 



keV counts s 1 10 6 counts s 1 arcmin 2 10 14 erg cm 2 s 1 dcg 2 kcV cm 2 s 1 sr 1 keV 



R12 0.11-0.284 0.025 6.85 4.66 0.55 

R45 0.44-1.21 0.02 4.76 31.3 0.83 

R67 0.73-2.01 0.02 4.91 26.9 0.42 



a Energy at 10% of the peak response. 
b Point Source Detection Limit. 
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Fig. 8.— The contribution of unresolved stellar 
sources to the RASS as a function of Galactic lat- 
itude. It is assumed that the point source removal 
limits were 0.025 counts s _1 in band R12 and 0.02 
counts s _1 in band R45 and band R67. 
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